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Abstract 



The Monte Carlo program, PYBBWH, is an implementation of the associated 
production of a charged Higgs and a W boson from bb fusion in a general Two- 
Higgs-Doublet model for both CP-conserving and CP-violating couplings. It is im- 
plemented as a external process to Pythia 6. The code can be downloaded from 
http: //www. isv.uu. se/thep/MC/pybbwh. 

1 Associated and W boson production 

The program, PYBBWH, is an implementation of the production of a charged Higgs boson, 
H^, in association with a W boson. The code was developed for the research presented 
in [1] where details on the theory and numerical results are presented. This program is 
written for a general Two-Higgs-Doublet model type II. The dominant production mode 
at tree level occurs is via bb fusion and at one-loop-level via gluon fusion. This program 
implements the leading order bb fusion part and the relevant Feynman diagrams are given 
in figure 1. By only including bb fusion this program is most suited for intermediate 
masses and large tan/3. 

In a general type II 2HDM the couplings relevant for this production can be specified 
via the Higgs mixing matrix Oji in the following way 

9h,h-w+ = 9HiH+w- = C>2i COS (3 - On sin /5 + i O^i , 
QHibh = Oii + i 03i sin p . 

In a real 2HDM were Hi = {hf^, H^, ^4°} the mixing matrix has the simple form 
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Figure 1: Feynman diagrams for H^W^ production via bb annihilation, bb — > H~^W~. 

which gives purely real couplings for h^, and imaginary couplings for A^. In a general 
2HDM there can be mixing between the CP-even and CP-odd Higgs states and the mixing 
matrix can have all elements non-zero. 

Using a formalism with only diagonal propagators for the Higgs bosons^, the differential 
cross-sections implemented in this program for the two processes are [2,3]: 
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where s and t are the ordinary Mandelstam variables of the hard process and 

X{x, y, z) ^ x"^ + y'^ + z'^ — 2{xy + yz + zx) , (5) 



^This formalism is valid as long as the Higgs bosons are well separated in mass. If they are so close 
that they overlap non-diagonal propagators are needed and this program can not be used. 
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with 9 being the polar angle in the 2^2 cms and 

1 

s — mi, + iniH^H 

Tit 13 

the propagators of the neutral Higgs bosons. 



2 _ A(s,m^,m|±)sin^6> 



2 Implementation 

The associated production program is implemented as an external process to Pythia 6. 
It uses the Les Houches generic user process interface for event generators [4] but it also 
uses Pythia specific routines. This means it can not be used directly together with other 
event generators, but after some minor changes to the code it should be possible. The code 
has been tested with with Pythia version 6.324 and the most recent version 6.413. Details 
on how to use an external process can be found in the Pythia 6 manual [5], section 9.9. 
In principle this is done by initiating Pythia with PYINIT( 'USER' , ' ' ,' ' ,OdO). The 
program does not work directly with Pythia 8 but it can be made to work via the Les 
Houches interface in Pythia 8. 

The widths of the and W bosons are included in the same way as in standard 
Pythia whenever possible, see details below. In other words the H"^ and W masses vary 
according to Breit-Wigncr distributions with mass dependent widths meaning that for each 
mass the decay widths arc recalculated based on the open decay channels. 

The program is setup for simulating LHC events as default, meaning incoming protons 
with 7 TeV energy. This can be changed by setting the parameters IDBMUP and EBMUP 
to incoming particle type and energy respectively. These parameters has to be set before 
PYINIT is called. 



2.1 Pythia specified couplings 

Pythia contains different SUSY simulations. The default in the PYBBWH program is to 
assume that one of these is used, so that the Higgs mixing angle a and tan /? are given in 
RMSS(18) and RMSS(5), respectively. The Higgs mixing matrix and the couplings used in 
the process generation is then calculated from a and tan p. Some checks are performed to 
see if the values of RMSS(18) and RMSS(5) correspond to the Higgs-fermion couphngs used 
in Pythia for Higgs decay. A warning is printed if an inconsistency is detected. 

2.2 User specified couplings 

If an external SUSY simulation or a general 2HDM type II is used, especially one with 
CP- violation, the default method can no longer be used. In this case the common block 
PYBBWH is used. It is defined as 
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INTEGER IPYBBWH 
DOUBLE PRECISION 0_M 
C0MM0N/PYBBWH/0_M(3,3) .IPYBBWH 

where 0_M(3,3) is the Higgs mixing matrix. IPYBBWH is a switch with the default vahie 
of 0, meaning Pythia specified couphngs. For IPYBBWH=1, 0_M(3,3) is used to calculate 
the couplings. Also in this case tan/3 is taken from RMSS(5). The masses and widths for 
all Higgs particles must also be set correctly in the PMAS array. 

Another difference between IPYBBWH=0 and IPYBBWH=1 is that for the former a varying 
width is used for the charged Higgs but for the later a fixed width, the one given by 
PMAS (37, 2), is used. This behavior is motivated since if the Higgs mixing matrix is 
specified in 0_M the couplings used in Pythia to calculate the Higgs decay are probably 
wrong. If ones wants to use a varying width in this case this can be done by substituting 
the function 

DOUBLE PRECISION FUNCTION HPWID(MHP) 

with a new function that gives the correctly varying width. The argument given to HPWID 
is defined as 

DOUBLE PRECISION MHP 

and is the mass at which the width is to be calculated. 

3 Download 

The code can be downloaded from http://www.isv.uu.se/thep/MC/pybbwh. On that 
web page there is also two example programs, one for Pythia specified couplings and one 
for User specified couplings. 

4 Final comments 

This code was developed for research published in [1]. The work was done in collaboration 
with Stefan Hesselbach and Johan Rathsman. If you use the code please cite [1] and this 
manual. 
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